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Abstract 

We study the thermal fluctuations of anisotropic order parameters (OP) in lay- 
ered superconductors. In particular, for copper oxides and a d-wave OP, we present 
some experimental consequences of fluctuations in the direction normal to the layers. 
It is shown that the c-axis penetration depth A c can have a " disorder-like" quadratic 
temperature dependence at low temperature. 

The fluctuations are analyzed in the framework of a Lawrence-Doniach model 
with an isotropic Fermi surface. Anisotropies pin the orientation of the OP to 
the crystallographic axes of the lattice. Then we study an extended t-J model 
that fits Fermi suface data of bilayers YBCO and BSCCO. This leads to a d- 
wave OP with two possible orientations and, including the thermal fluctuations, 
yields the announced temperature dependence of A c . Furthermore a reservoir layer 
is introduced. It implies a finite density of states at the Fermi energy which is 
successfully compared to conductance and specific heat measurements. 
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1 Introduction 



Despite an enormous amount of theoretical and experimental research devoted to 
High Temperature Superconductors (HTSC) in the past decade, many important 
issues still remain unsettled. In particular there is as yet no explanation for the 
pairing mechanism at work in these materials. Early on Anderson suggested [ I J 
that spin fluctuations might play a critical role for the superconductivity of copper 
oxides. These fluctuations have been analyzed for instance in the framework of the 
t-J model [0 and a natural candidate to represent the superconducting state in that 
model is a d-wave order parameter (OP). 

Many low temperature experiments seem to suggest that a d-wave symmetry is 
a strong possibility : the variation of the penetration depth with T |, | and H 
[|], Knight Shift and \jT{T data |J, angle resolved photoemission spectra f7|, the 
variation of Cy with T and H ||. Yet these experiments do not probe directly 
the OP so that one cannot rule out other symmetries or mixtures of symmetries 
JJ. Precise measurements of the single particle density of state (DOS) near the 
Fermi energy in the superconducting state ought to give some clue about the pairing 
symmetry. Tunneling experiments however have been fairly inconclusive thus far, 
leading to s-wave, d-wave and, in the case of YBCO, to the conductance of a 
system with finite DOS at the Fermi surface [K], p| . Materials issues have made it 
difficult to sort out genuine features from experimental artefacts. The observation 
of spin fluctuations well below T c seen by inelastic neutron scattering and a linear 
in T dependence of Cy once the phonon contribution is subtracted ||T2[ have been 
interpretated as evidence for gapless superconductivity. 

Despite the fact that a clear-cut answer still remains elusive for high T c materials, 
an interesting by-product of these studies has been a renewed interest in coherence 
effects in d-wave superconductors. Sigrist and Rice indeed suggested that d-wave su- 
perconducting currents would display a strong spatial anisotropy at a boundary |13| . 
This anisotropy should show up in Josephson experiment setups [fL4] . Ambiguous 
results have been obtained, although Tsuei-Kirtley's experiment strongly supports a 
scenario where one has nodes for the OP (consistent with either d-wave or extended 
s-wave). Furthermore, the layer structure caracterizing C11O2 compounds should 
affect tunneling and Josephson currents as it does for vortex (thermo)-dynamics 
below T c and for superconducting fluctuations above T c , where a Lawrence-Doniach 
model is used [jn|. These reasons led us to investigate the issues of layering and 
of anisotropies in d-wave superconductors. Our results show that these effects may 
play a significant role in the interpretation of tunneling and penetration depth types 
of experiments. 

Our paper is organized as follows : in section 2.1 we present the Landau-Ginzburg 
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free energy of a single layer of a d-wave superconductor with an isotropic Fermi 
surface by solving the Gorkov equations. We calculate the energy cost of the fluc- 
tuations in the orientation of the nodal lines, denoted a, of an OP with the d x 2_ y 2 
symmetry. These fluctuations can be seen as admixing d x 2_ y 2 and d xy terms and 
appear in gradient terms of the free energy. We also show that the d-wave symmetry 
results in anisotropics in gradients of a and of the phase of the OP. These do not 
modify the critical behaviour near T c but have an influence in the London regime at 
low T. Then, introducing anisotropies in the Fermi surface, we see that this tends 
to orient the nodal lines at 45° with respect to crystallographic axes. We next con- 
sider the case of weakly coupled d-wave superconducting sheets and we derive the 
corresponding Lawrence-Doniach free energy (see section 2.2). The energy cost of 
the fluctuations of a is once again computed and gives, in particular, the Josephson 
current perpendicular to the layers. We find that the fluctuations in the c-direction 
are limited to a coherence length. 

In section 3 we consider a t-J model that reproduces a reasonable band structure 
for the bilayer compounds YBCO and also BSCCO. It is found that the nodal lines 
can now take on two positions This seems in agreement with recent experiments 
on BSCCO [[17). We add a hopping energy t n that couples two bilayers via a charge 
reservoir; depending on the strength of t n one may go from a more 3D structure 
pertaining to YBCO to a more 2D structure relevant for BSCCO. We show how 
these modify the DOS which is determined analytically near the Fermi surface. It 
starts linearily, as in case with pure d x 2_ y 2 OP, with an offset that vanishes with t n . 
This offset is observed in specific heat measurements (l8j that are in quantitative 



agreement with our result and also in conductance measurements [|H], |TT[ . 

The existence of several orientations of the nodal lines enhances the effect of 
the thermal fluctuations pointed out above: they give a mechanism that permits 
the OP to change its orientation from one layer to another. This has experimental 
implications that we discuss in section 3.3, in particular in the case of c-direction 
penetration depth measurements. Fluctuations in the orientation of the OP in the 
c-direction can change the T dependence of A c from a linear to a quadratic behaviour 
at low T, even in pure samples where it is linear in the a- and b-directions. This is 



consistent with experiments on YBCO 
Section 4 is a brief conclusion. 
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2 Fluctuations of d-wave order parameters : phe- 
nomenology 



The model under consideration is the one introduced by Lawrence and Doniach 
(LD) |15|] for weakly coupled layered superconductors. Here we shall first study 
the motion of electrons within a layer and then consider the Josephson tunnel that 
couples adjacent layers. 

We introduce the "nodal angle" which define the orientation of the nodes of the 
d x 2_ y 2 OP within a layer. It fluctuates both in the direction parallel and perpendic- 
ular to the layers and we determine the length scales of these fluctuations. 



2.1 Free energy density in a layer 

We consider a d-wave OP A n (k, R) where R is the center-of-mass coordinate of a 
Cooper pair, n the layer's number and k the relative coordinate in the Fourier space, 
k = (k x , k y ) = (|k| cos#, |k| sin#). We choose for the k-dependence of A the form 

X V 

= cos(26') which has nodal lines, k x = ±k y . If a n (R) denotes the relative 
angle of these lines with respect to fixed axes one has 

A n (k, R) = A (T) cos[2(fl - « n (R))]e^( R ) (1) 

This can be expanded as 

A n (k, R) = A?(R)V>!(0) + A£(R)V> 2 (0) (2) 

where the A" are two related gap functions 

/ A?(R) = A cos(2«) e^ 

1 A"(R) = A sin(2a) 1 ' 



and the ip™ have the form of the d-wave pairing called d x i_ y i and d xy 

(4) 




= cos (29) =k 2 x -k 2 y 
= sin(2#) = 2k x k y 

This is compatible with the factorized pairing potential V(k— k') —VJ2 ip i (li)ip i (k'). 
The form (§) for the OP can be generalised to other symmetries or mixture of 
symmetries. A relevant example is the case with a s + d OP (e.g. when one has 
orthorhombicity) which is given by @) plus a term A,^ip^ having the desired s-wave 
symmetry. We mention later the resulting corrections to the free energy. 
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Then we shall write the GL equations for A n where one integrates over k. We 
follow a method analogous to the one used previously for mixed OP |21| except 
that we introduce the nodal angle a n (R) and we consider a different set of A*}. 
The derivation of these equations may be done in two steps. The first concerns any 
OP and generalizes the s-wave pairing equations. Since this has been shown several 
times (see e.g. |]22|| ) we only sketch the calculation in Appendix A. The second step 
specializes to the d-wave pairing considered above ([I]). The whole procedure gives 
two coupled gap equations for the A" in which, in the limit T — > T c , we retain only 
the following terms A, V 2 A, V 4 A and A 3 . We get 

A? = a<V?>A?-A< V ,4 > | A „|2 A n_3 b <^2 > | A n|2 A n 

c c 

+V v 2 F (< tfkl > V 2 .+ < fift > V 2 ) A? 

+c v% (< ^ > Vt+ < ^XK > V: + 6 < i&lkl > V 2 V 2 ) A? 

+c v% (4 < ^ 2 k 3 x k y > V 3 V y + 4 < ^ 2 k x ky > V^V 3 ) A^ (5) 

In (||), < ... > denotes an integration over 6 1 and the various letters a, b, denote 
the integration and summation / d\k\ J2u ■■■■ The Fermi surface being circular the 
energy only depends on |k| and the integration on 9 affects only the i\> i and the Fermi 
velocity [i.e. (k x , k y ) = mvF{cos9, sin^)]. One gets a similar equation, upon making 
the substitution A 1 <-> A 2 (and ip 1 <-h> ijj 2 ), which has the same coefficients a, b, rj and 
c. In these equations the d-wave pairing manifests itself in the < ipl >-dependence 
of the coefficients (while a non zero value of a n introduced two gaps A"). 

It appears that the contribution of the fourth order in V is anisotropic (see the 
c-terms in ((D). In the case of s-wave pairing these terms can be combined to give 
a single isotropic contribution to the free energy, as it should. 

Finally the GL free energy density of the nth layer is given by 



fn — fo + fani (6) 



with 



k = -4(1-^)(|A?| 2 + |A?| 2 ) + ^ (|A?| 2 + |A™| 2 ) 2 



+^(^) 2 (|V R A 1 | 2 + |V R A 2 | 2 ) 

J- c 

and the anisotropic part 

fam = E (^) 4 [ C t (|V 2 Alf + |V 2 A™| 2 ) + C[ | V.V.Alf 

i=l,2 ±c 
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+- (^) 4 1 v,a?v^t + v,Arv;A; - v^v^r - v^rv^ 

+(Aj «-> A 2 ) | 

where we used @ and the circular model to calculate the integrations over 9 which 
we then incorporated in the coefficients A,B, .... 

This free energy density has been obtained keeping only the leading term in the 
small- 1 k | expansion which yields the from ([!]) for the OP and an isotropic Fermi 
surface (FS). Within this restriction we got an expression that does not depend 
on a n (R) at the mean- field level. In other words, the free energy of one layer is 
minimised by any uniform a n . 

Let us now consider the anisotropy effects. They are two sources of anisotropy, 
one in the OP, the other in the Fermi energy. The effect of the first one is illustrated 
by the presence of f an i which implies that the fluctuations of a n (R) (and also of 
n (R)) are anisotropic in a layer. The fact that this anisotropy appears to order 
Vr 4 A 2 is due to the particular form of the OP (another form yields another order, 
e.g. an additional s-wave OP gives anisotropic terms to order Vr 2 A 2 0). 

One can verify that the anisotropic terms f an i contain gradients of a and of 
cf) which are coupled (this coupling being absent at lower orders). These terms are 
negligible in the high temperature region (near T c ) but come into play at lower T. In 
the London regime, it could be interesting to study the implication of this coupling 
on the vortex state. As one replaces gradients by V — ^ A the term to the quadratic 
order takes on the form 

A*£{(^) 2 [4(V R a B ) 2 + (A + V R n ) 2 ] (7) 

In addition to the kinetic energy of the current one gets a rigidity term in the 
orientation of the OP, a n . Thus the fluctuations of a n are relevant as long as a n is 
not fixed relatively to the crystallographic axes. 

Indeed the introduction of an anisotropy in the FS due to a crystal structure can 
pin in the orientation a n and also generate anisotropic fluctuations of a n . To show 
this on an example, we consider an elliptic FS, i.e. fc 2 + /c 2 + e/c 2 . To first order in e 
the anisotropic terms now come in to order V 2 A 2 while to order e 2 the terms |A™| 2 
and (Agl 2 have different coefficients. This means that the term to order |A| 2 now 
depends on a n which will be fixed by the minimisation of the free energy. In this 
case there are still anisotropic terms in the Free energy (in gradients of the phase of 
the OP) which are relevant in the London regime. 
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2.2 Free energy density for layered superconductors 

Following \TE\ we now add a transverse component to the electron energy of the 
normal state 2 

e (k) = k -l±^L + K cos ( kz d) (8) 

where d is the interlayer spacing. We performed an expansion of the gap equations 
in the limits T — T c << 1 and K « T c (see Appendix A). The free energy has the 
general form AF = £„ J dR(f n + f n>n +i) with /„ given by (§) and 



/n,n+i = ^(lA^-Ar^ + IA^-A^ 1 ! 2 ; 



+12 C7 ^ (^) 2 (|V R (A? - A™ +1 )| 2 + |V R (A« - A« +1 )| 2 ) 

+^2 C jf £ (isK-Ar 1 ! 2 

- 8|A™ - A™ +2 | 2 + 3|A^ - A™ +3 | 2 - ^|A™ - A^ +4 | 2 ) (9) 

Thus in this LD model thermal fluctuations give rise to variations of the value of 
the nodal angle a n . One sees in (^|) and ([]) that a n fluctuates, both within a layer 
and from one layer to another, around a constant but arbitrary value (in f the first 
two terms do not depend on a n ). The energy cost of these fluctuations is given by 

5f = A 2 E {4 (^) 2 (V R a n ) 2 + 2 ^ (l-cos[2(a n -o n+1 )] cos(0 n -0 n+1 ))} (10) 

which was determined from the gradient terms of /„ and j n . n +i with the simplifying 
assumption that Vr0„ = 0. We wrote only the dominant terms, that is, the first line 
in @ and the part fo of /„ but one must keep in mind the presence of anisotropic 
terms, from ([It]) one gets a current density between two adjacent layers 

j n>n +i ~ cos[2(a n - a„ + i)] sin(0 n - (p n+1 ) (11) 

Although this is a current that flows perpendicularly to the layers, it has the same 
behaviour with a n — a n+ i as the current across a tunnel junction in a layer |23| . 
Note that in a model where a n — a n+ i can become arbitrarily large (i.e. when a n is 
not pinned by anisotropies in the FS), the fluctuations can supress the tunneling in 
the c-direction (j n ,n+i vanishes on average). 
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Now to study the energetic balance of the fluctuations of a n one has to compare 
the two characteristic lengths, £y and that come from 5f and a third one, £ an ,, 
from f ani : 



£ ± 2 K 2 d 2 



^ 2 (/^ 



(12) 



and 



where AC/E 2 is a numerical factor. The LD model is valid within a regime of 
temperature given by the limits 

1 - £ « 1 and a ~ (1 - ^r V2 « d (14) 

J-c *-c 

There, the intra-layer length £ ani is of order of £y, while the inter-layer length, 
is much smaller yielding the fluctuations of the OP in the c-direction inexpensive. 
We study some experimental implications of the fluctuations in the c-direction in 
section 3. 

The above analysis does not fix the value of a n in the absence of anisotropics of 
the FS linked to the crystallographic structure of these materials. To remedy this 
one has to consider a microscopic model that specifies a reasonable band structure 
(a step in that direction was taken in section 2.1). We choose to start from the t-J 
model since it is a natural candidate for a d-wave superconducting phase Hj. 



3 The t-J model for bilayers coupled via a reser- 
voir 

In the usual t-J model one considers the pairing of nearest-neighbor fermions on a 
square lattice with an exchange energy Jn and a quasiparticle hopping energy t\\. 



This leads 0, ^6] to a normal state with a Fermi surface (FS) of square shape at 
half filling and favors a variational d-wave superconducting state (in fact the lowest 
energy state corresponds to an s+id OP at half filling and to a d-wave OP away from 
half filling) There the OP has the d x i_ y i symmetry of flU) with a nodal angle a = 0, 
( nodes lie at 45° with respect to the crystallographic axes ). Recently it was shown 
[ |j~6|j that additional terms were needed to fit experimental data on BSCCO: Kuboki 
and Lee used a bilayer structure with an exchange energy between two layers, J±, 



and a hopping term, t±. In agreement with ARPES experiments one gets a FS 



with two bands [IB|. Furthermore the OP has two nodal angles at a = ±«o which 
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may have been observed in BSCCO jI7| («o ~ 10°). Here we shall go further by 
adding terms that describe the FS and the density of states (DOS) of YBCO. We 
include a pairing along the diagonals in a layer, J^, and a hopping term, t n , that 
couples two bilayers via a charge reservoir, see figure 1. 



3.1 Fermi surface and gap 

The slave-boson mean-field Hamiltonian of the t-J model describes electrons with 
forbidden doubled occupancy on a site : the creation of an electron on site i in 
the layer s consists of the anihilation of a hole operator having the charge of the 
electron 6| and the creation of a fermion operator having the spin a of the electron 
fi,a ( s = for the bilayer and s = 3 for the reservoir). To simplify we shall set 
b*l + bj = 5, i.e. the number of empty sites gives the doping rate. 
The hamiltonian is given by 

n = -EEt 5 % f&fja + Jv - \n inj ) ] (15) 

with the constraint fujfio = 1 — 5 and S? = fia&apfjp, the Uj and are shown 
in figure 1. Within mean-field approximation we introduce various OP, in particular 
we restrict our analysis to a d-wave coupling in the layers and to an s-wave coupling 
between layers : 



J\\ E °- < flafj-a >= +(-) A|| with j=i± x(±y) and r = 1, 2 
a 

J dJ2 ° < flaf 3 -a >= +(-) A d with j = % ± x ± y(T% =F y) and r = 1, 2 
a 

J±J2 ° < fiaff-a >= A ± with j = i + z 



(16) 



in the particle-particle channel and 

E < fia fh >= X\\ with j = i ± x or ±y 
a 

E < /icr fja >= Xd with j = i ± (z ± y) 
a 

E < /«t //a >= X± with j = i + z 



(7 



E < flafjcr >= Xn with j = i + 



IjCT 

a 



(17) 
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in the particle- hole channel ( here too r = 1 , 2 ) . 

As usual H, |26| we write the Hamiltonian terms depending on the OP in mo- 
mentum space 



3N Ajj 



A* A 1 



H A = £F+(k).M(k).F(k) + ^(2^ + 2 

with F(k) = (/£ T , /l k)i , £ T , /! k4 , / k 3 T , /! k J, 



M 



J 

\ 



' ii 




/ £l Ai t± A12 tn 

-6 A 12 -t ± -/ 

-6 o -t n 
6 o 
-6 / 



V 

where one easily completes the matrix knowing that M t - 
gaps written below (taking the lattice parameters a = b 



M* and the energies and 
c=l) 



6 
£ 3 



— e || (cos k x + cos k y ) + 2ed cos k x cos k y — jj, 
— e n (cos fca; + COS fcj,) — jtt 



-e±e 



\k z d 



(19) 



with ta = St® + <^aXa f° r the various energies and /i is the chemical potential which 
depends on 5. With these expressions one can reproduce the FS of YBCO while 
the FS of BSC CO is well represented if one neglects the next-nearest interactions, 
e d = 0, and the hopping between bilayers, e n = 5t n = 0. The FS of figure 2 is 
obtained upon adjusting the various parameters t a and J a in order to reproduce 
experimental data on optimally doped YBCO (i.e for 5 ~ 0.25 ) P5] , that is (in 
eV), t\\ ~ 0.4, t d ~ 0.1, t± ~ 0.05 and t n ~ 0.0045 . (We determine these values 
from t n = T 2 (a/d') 2 , and t± = T 2 (a/d) 2 , using T = 1/5, d'/a = 3,d/a = 1 while td 
comes from [p7[| .) As it has been pointed out in case of BSC CO [0 the hopping 
term between layers is taken in the form t± ~ (cosk x — cosk y ) 2 . 
The gaps in M are given by 

Ai = A||(cosfc x — cosky) + 2A d sm k x sin k y 



A12 = 



A ± e ifczd 



(20) 
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so that the pairing along the lattice axes is of d x 2_ y 2 type and, along the diagonals, 
of d xy type. 

One can diagonalise the Hamiltonian (|i8|) which gives the quasiparticle energies 
E h2 = ±[(ei-t±) 2 + (A 1 -A 12 ) 2 ] 1/2 

£ 5 , 6 = ± [(e 1 + t ± )2 + (A 1 + A 12 ) 2 + e 3 2 + 4t 2 n -5 1 / 2 ] 1 / 2 (21) 

with 

S = [(6+^) 2 +(A 1 +A 12 ) 2 -4 2 ] 2 +8t 2 [(6+^) 2 +(A 1 +A 12 ) 2 +2e 3 fe+ei+^)] (22) 
This yields the following free energy 

Q/V Am A 2 A 2 

F A = -2T££ln[2cosh(£ r (k)/2T)] + ^(2-1 + 2-* + -±) (23) 

r , 8 «V|| Jj_ 

One then minimises (^) which gives the various OP. Here we restrict ourself to 
the above forms, see fl2"0|), for the OP, i.e. the degrees of freedom are An, A^ and 
A^ which are real. There are three critical temperatures repectively proportional 
to J||, Jd and Jj_; with the values that give figure 2, J\\ dominates. In analogy with 
our previous OP, see (|l|), we set 

f A,, = A cos(2a) 
\ A d = A sin(2a) 

Clearly the value a (mod 7r/2) = (i.e. nodal lines at 45° with respect to the lattice 
axes) will minimise the free energy since Jn » Jd (this can be verified e.g. upon 
expanding the free energy in the small- A limit). The additional coupling Jd creates 
a metastable state in which the OP has nodal lines along the lattice axes, separated 
from the lowest energy state by an energy barrier of order Jn. 

From (pip one easily verifies that the mixing of Ai and Ai 2 leads to two nodal 
angles obtained by setting Ai ± A i2 = and denoted by a = ±a . Considering 
the effects of the thermal fluctuations on the layered system studied in the previous 
section one gets an OP whose orientations fluctuate from one layer to another be- 
tween the values a = ±a - This will have dramatic consequences discussed in the 
following. 
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3.2 Density of States and conductance 

In order to compute the DOS for this system one must have the eigenvectors of the 
matrix M. Indeed, denoting by U the unitary matrix which diagonalizes M, the 
DOS in the s-th layer is given by 

p s {E) =2j2\U s ,i\ 2 S(E-E l (k)) (25) 

k,Z 

where I indexes the eigenvalues and E is the energy measured relatively to the Fermi 
energy In the absence of inter-bilayer coupling the Ei(k) yield "nodal" contribu- 
tions for the DOS in layers 1 and 2 and "normal" contributions for layer 3. A nodal 
contribution is obtained whenever both the £ and the A terms vanish in the quasi- 
particle spectrum; in that case the DOS varies linearly with E for small E. If only 
a £ term vanishes, layer 3 which has A = yields a finite contribution for vanishing 
E. Including the inter-bilayer coupling will admix these two types of contributions 
such that by summing over k in any given layer one may obtain a linear combina- 



tion of "nodal" and "normal" DOS [28]. ( As noted above the nodal and normal 
contributions arise from distinct ranges of integration in k space ). 

With our model it is straightforward to compute IA and we only quote the result 
for the DOS in layer 1 : 

Pl (E) = aE + b (26) 
where b is proportional to the normal state DOS in layer 3, pjy, 

b = PN^r (27) 

and 

„&) = 2 6 + fx-(teH-^ + (A 1+ A^)^ (28) 

Ai + A12 

with k such that £s(k ) = 0. The offset at the Fermi energy is also present in a 



DOS determined from an s-wave OP p8| while the linear in E dependence is due to 
the d-wave symmetry. 

In single-crystal YBCO, the DOS has been determined from the specific heat 
[ ITH ] which is dominated by a T-term at low temperature whose prefactor yields 
Pi(Ef)/pn ~ 1/5. This is in close agreement with what we get for (^) 2 using 
t n = 4.5 10 -3 eV and A YB co = 15meV P|. On the other hand we get for BSC CO 



Pi(E F )/p N ~ 1(T 8 using t n = 4 10" 6 eV and A BSC co = 20meV. 

Furthermore the DOS (EBf) yields a tunneling spectrum along the c axis shown 



in figure 3 at zero temperature from 

d r eV 



G{ y) ~ -ITr I dEp (E) Pl (E - eV) (29) 
aV Jo 
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where po is the DOS of a conventional superconductor. The predicted spectrum 
closely ressembles what is found in experiments on the tunneling between YBCO 
and Pb [|T0 , |TT| . Above the critical temperature of Pb the finite offset at zero voltage 



that is measured is given in our analysis by the hopping energy t n (see (0)). Below 
this critical temperature, the offset induces the divergence at A Pb = 1.38 (see figure 
3). Furthermore there is possibly a linear regime in these experimental data coming 
from the d-wave symmetry of the OP of YBCO. Finaly one expects a slope change 
at eV = A Pb + Aybco- ( m figure 3 we use A Y bco — 20meV which is lower than 



above since we now take the value given by conductance measurements [|T0| , |TT 



On the other hand no offset seems to be observed in experiments on BSCCO [29]. 



3.3 density of states and penetration depth in the c-direction 

For a 2D d-wave superconductor one finds that the DOS depends linearly on E for 
small E (see e.g. pO] for a complete result). Indeed, 



p(E) oc J J kdkde 5(E - Jil + Al) (30) 

where Ak = Aq(T) cos(2#). We now compute the contribution to p near 9 = 
\ (mod |). Setting 9 = f - 9', x = £, y = 2A 9', r = ^x' 2 + y' 2 and ip = tau -1 (y/o;) 
we obtain 

p(E) oc J J rdrdi) 8{E - r) (31) 

Thus the linear behaviour is due to the nodal structure of Ak which lead us to address 
the question of the effect of the thermal fluctuations on the DOS and consequently 
on the penetration depth. 

As discussed in the previous section, for bilayers, the contribution at 9 — f 
(mod |) evolves into two modes at 9 = | ± a (mod |). At T = the inter-bilayer 
hopping term t n couples the bilayers and selects the +a Q or — a solution to form a 
coherent state. For T > t n however the bilayers are decoupled. The system in effect 
behaves in a 2D fashion. As a result, the contribution to the DOS along c coming 
from a particular orientation of the nodal line - say 9 = f — «o - is given by 



p c (E, z) oc I dip rdr 5(E - r) (32) 

Je{z) 



with e(z) = or A ao so that 

p c (E, z) oc E for E > e 



p c (E,z) = for E<e ^ 
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The average DOS is thus given by p c (E) = < p c (E,z) > e where < .. > e denotes 
averaging over the e distribution function. Since the nodal lines may take on two 
positions the system is equivalent to an Ising problem. The energy barrier between 
the two equilibrium positions in a bilayer is of the order of A ao, (and can be 
estimated to ~ 40K using realistic values for A ~ 25meV and for a ~ 10°). 
The energy barrier in the vertical direction is of the order of t n . We note that for 
T > A a contributions to the DOS are not dominated by the nodal lines anymore, 
and gapped parts in k space must be included. In that limit the Ising analogy ceases 
to be valid. In the Ising regime (T < Ao«o) however we may write 

p(£(z)) = (34) 

where the variance a oc ^ (the susceptibility divided by the linear size along the c 
axis). Thus 

p c (E) = -== / e^de (35) 

V27TCX 2 Jo 

If the susceptibility is large enough, as is the case for a quasi 2D system and/or if L 
is small enough as is the case for thin films, one may have a situation where E oc T 
and a » T so that 

Pc{E) oc E 2 (36) 

This situation may correspond to high T c films and to the bismuth, mercury, and 
thallium compounds which behave as fairly 2D systems. On the other hand for large 
L and/or small x as is the case for single crystals and/or more isotropic compounds 
such as YBCO one recovers 

Pc(E) oc E (37) 

Inasmuch as in the thermal fluctuation dominated regime the penetration depth is 
given by 



K 2 



where f(E) = exp ^ E ^ +1 and Ptv(0) is the normal state DOS at the fermi energy, one 
sees that thermal fluctuations may genuinely yield 

AA; 2 oc T 2 (39) 

even in the absence of disorder 



Experimentally it seems indeed that for the 2D-like materials a T 2 law is ob- 
tained for the c-axis penetration depth. For YBCO the situation is still not settled. 
Mao et al. recently reported findings of a linear T dependence for a YBCO single 
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crystal || whereas the UBC group seems to find a more quadratic dependence |H5| . 
Moreover experiments on YBCO powders also see a T 2 dependence that is shown 
to be intrinsic |20[| . As was pointed out in |2(| the measured Ay combines \ a b and 
A c which is much larger. So even if the present model would give a linear behaviour 
for AA a &, one could find AAy oc T 2 in case of a quadratic AA C . 



4 conclusion 

In this paper we consider superconductors consisting of layers with a d-wave OP 
that varies in momentum space according to 

A k = A (T)cos[2(#-a)] (40) 



(tan# = ^-) where a gives the nodal lines orientation. As for the phase of the OP, 
a may fluctuate in space. We consider an a that is homogeneous in the layers but 
can thermally fluctuate in the c-direction. We show how the value of a is fixed 
with respect to crystallographic axes when a reasonable band structure is taken into 
account ( this leads to anisotropics of the Fermi surface ). We use an extended t-J 
Hamiltonian with Cu02 bilayers separated by a charge reservoir, yielding a FS for 
YBCO in fair agreement with that mapped out by ARPES measurements. This 
leads to various possible pairings of the type of ( f40"D with different a. Using the 
key idea of this paper of combining the inexpensive spatial fluctuations of a with 
the fact that the OP can choose between several orientations, we show that the 
thermal fluctuations of a change the T dependence of AA C from linear to quadratic. 
Our model predicts this quadratic behaviour for thin films and quasi-2D systems 
where it is usually attributed to the disorder (see |32| for a review). We also show 



how the presence of the reservoir layers is revealed by conductance and specific heat 
measurements. 

Here we did not consider the orthorhombicity that distinguishes YBCO, from 
e.g. BSCCO, that is, asymmetric d-wave amplitudes in the a and b directions. 
Then one gets a mixture of d x 2_ y 2 and s x 2 +y 2 (called extended s-wave) symmetries. 
If the s-wave component is small one still has node lines. In twined samples there 
will be domains of different orientation of these lines, an issue that we leave for a 
later publication. 

Finaly let us consider the case of a pure extended s-wave OP sometimes invoked 
for BSCCO. It has nodes and one can easily show that they yields a linear behaviour 
of the DOS near the Fermi energy. Since there are thermal fluctuations in the 
orientation of the OP as well (see section 2), one can conclude that the above results 
for the penetration depth persist. 
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Appendix 

In the following we give a summary of the derivation of the LG-equations for a 
general OP in the vicinity of T c . As in section 2 we first derive the equations in a 
layer and then consider the LD-model. The OP is defined by 

ZKry)=TV>-r')£F+(r,r') (41) 

where V is the attractive two-body interaction of the weak-coupling theory and 
uo = irT(2n + 1). We start from the Gor'kov equation of motion for the normal and 
anomalous Green functions (in the absence of a magnetic field) 

/ (iw - k') + E q 4(k, q)F+(q, k') = 5 2 (k - k') 

1 (^ + e k )i^(k,k / )+E q ^(k,q)G c ,(q,k') = ^ 

where we omit the spin variable to simplify. The energy is described by £(k) = 
e(k) — fi with e(k) = {k 2 + ky)/(2m). As usual one introduces the Green function 
of the electron in the normal state 

^(k,k') = ^^ (43) 

and one expands in power of A. One gets 



A*(k u k 2 )=TY,H <^( k i + k 2 - k - k')V(k - kx) 
u k,k' 



A*(k,k') 



(iuj - &)(-«<; - Ck') 



A*(k, q )A( q , q >)A*(q,k>) 1 
^ {-iu - f k )(iw - i^{-iuj - f q 0(iw - J 

We shall write Z\(ki,k 2 ) = A(ki — k 2 ,ki + k 2 ) which is the notation of section 2. 
Then we make the substitution in various quantities 

q 2 k Ikl 

6c+ q = £k + |u|q.k + — with k = — , \v\ = — (45) 

Ann LK Tfl 
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and we expand in powers of q (we set = iuj — £k) 




1 

iu - £ kH 



In the following we shall keep only the dominant terms without loss of generality. 
This becomes valid when uj 2 « /j 2 which leads to substitute factors Vp = (jj/2m) 1 ^ 2 
for \v\. Furthermore one substitutes the integration m J d^J dd for J2k t° which 
the terms with an even power of k = (cos 9, sin 9) won't contribute. Thus we only 
keep in fl46|) three terms, the order in q and the first term written above at the 
order 2 and at the order 4. 

Then the particular forms of the OP and the pairing potential, see ([!]), (|2|) and 
below (|J), are to be introduced. Upon taking the Fourier transform with respect to 
the center-of-mass coordinate q, which gives gradients Vr of A, one gets ((5|). 

In case of layered systems the gap equations can be derived as above, with the 
energy In the expansion of (|4"3| ) we keep at each order in (q, q z ) the dominant 
terms in the limits \u\ ~ T c << \i ~ vp and K << T c at finite K 2 /(T c v F ), that is, 



I 1 L \v\ 2 , * 2 K 2 




QvlK 2 



q.k) 2 + — [l-cos(dfe)] 



r)l vq-k) 4 + ^^(q-k) 2 [i-cos(^)] 

3K 4 15 1 1 

+ 16£> 2 _ 13cos(<ig 2 ) + 8cos(2dg 2 ) - 3cos(3dg 2 ) + -cos(Adq z )\ j (47) 

As before we also anticipated upon supressing terms that disapear after an integra- 
tion over (|k|k, k z ). One substitutes (f37D in ( fH|) to get the gap equations and finaly 
the free energy densities (0) and (P). 
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Figure caption 



Figure 1 : 

The hopping and exchange energies on a bilayer structure (1-2) separated by 
charge reservoir (3). The summation in the Hamiltonian (15) is restricted to nearest 
and next nearest neighbors witin a layer (right figure) and to nearest neighbors in 
adjacent layers (left figure). 

Figure 2 : 

The Fermi surface of YBCO7 for a single bilayer in the TXY quadrant. The 
solid (resp. dashed) line corresponds to a dispersion relation of the form £1+ (resp. 
- ) tj_ (see eqn (18)). The degeneracy of the two branches at 2 is lifted by t n . 

Figure 3 : 

Zero temperature tunneling conductance along the c direction for a YBCO/Pb 
junction. The DOS for YBCO near the Fermi energy is given by eqn (21). The 
conductance diverges at the Pb gap taken to be 1.38 meV. 
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